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1. Introduction 

Lattice calculations are becoming increasingly important for particle physics phenomenology. 
They address and quantify the "hadronic uncertainties" which still afflict many quantities that con- 
strain the validity of the Standard Model [1-6]. Almost every talk on lattice QCD delivered to a 
more general audience during the past 10-15 years contained the phrase that "lattice calculations 
are performed at unphysical quark masses." What we usually mean by this statement is that for any 
given discretisation, there is a priori no way of knowing which values of the bare quark masses 
correspond to those of the physical quarks. In most, if not all, cases it turns out that the physical 
light quark masses lie outside the regime which is directly accessible using the currently available 
algorithms and machines. Likewise, the physical values of the heavy quarks are dangerously close 
to, if not above, the affordable cutoff scale. 

Chiral extrapolations are thus required in order to make contact with the physical light quark 
masses. Chiral Perturbation Theory (ChPT) provides theoretical constraints on the quark mass 
dependence of observables, based on the underlying dynamics associated with chiral symmetry 
breaking. The following quotation from the FLAG report [1] serves as a reminder that improving 
the control over the chiral behaviour is mandatory in order to make further progress: 

"Although [light quark masses] are decreasing very significantly with time [. . .] it remains 
true that [the chiral] extrapolation is one of the most significant sources of systematic error." 

Some collaborations have produced lattice data at or even below the physical pion mass [7, 8] 
which may render extrapolations guided by ChPT soon obsolete. Instead, one can resort to some 
sort of analytic ansatz to interpolate lattice data to the physical point. However, even if ChPT were 
not required to perform chiral extrapolations, a comparison of the quark mass dependence deter- 
mined on the lattice with the predictions of an effective theory would provide useful information, 
since it allows for the determination of the low-energy constants (LECs) which parameterise ChPT. 
Furthermore, as simulation algorithms still show a significant increase in computational cost when 
the light quark masses are tuned to their physical values, it is still difficult to disentangle the chiral 
behaviour from systematic effects arising from finite volume and/or coarse lattice spacings. 

The vast majority of lattice estimates for phenomenologically relevant quantities is still dom- 
inated by systematic errors. In this review I try to investigate whether chiral extrapolations are 
sufficiently well controlled in order to rightfully claim the accuracy which is quoted in the recent 
compilations. Here I shall focus on three different types of observables: Lattice estimates of the 
light quark masses are discussed in the next section. In section 3 I will study the systematics of 
chiral fits applied to meson decay constants. Section 4 contains a discussion of the chiral behaviour 
of the axial charge of the nucleon. Summary and conclusions are provided in section 5. For the 
purpose of this review all lattice results are taken at face value. Issues such as "rooting", induced 
non-localities, or the freezing of topology will not be discussed here. 

2. Light quark masses 

Quark masses are fundamental parameters of the Standard Model whose values determine 
many important quantities in particle phenomenology. The recent compilation of lattice results for 
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Figure 1: Results for the strange quark mass in the MS-scheme at 2GeV, obtained in lattice QCD with 
Nf = 2 + 1 and Nf = 2 flavours of dynamical quarks [1]. Green points represent lattice results which are free 
of any red tags according to the FLAG criteria. Blue circles denote the results from sum rule calculations. 
The grey band and the vertical dotted lines denote the global estimate for the Nf = 2 + 1 and two-flavour 
theory, respectively. The PDG estimate is shown at the bottom. 



the light quark masses and their conversion into "global" averages in the FLAG report [1] is based 
on a set of "quality criteria". Using a simple colour code, they are meant to assess the quality of a 
given calculation regarding a number of different systematic effects. Obviously, these criteria must 
be adjusted over time, in order to reflect the true state-of-the-art. In the current FLAG review, a 
green star (+) is awarded if the systematic error is "convincingly shown to be under control". An 
amber ball ( ) signifies that a "reasonable attempt" at estimating a particular systematic error has 
been made. Finally, a red box (■) indicates that no attempt was undertaken to quantify a systematic 
effect. To set the scene for the discussions to follow, the criteria for the colour code referring to 
chiral extrapolations and the related finite-volume effects are repeated here: 

Chiral extrapolation: Finite- volume effects: 

★ m™' m < 250MeV ★ m™'"L > 4 or at least 3 volumes 

250MeV < m™ in < 400 MeV • m™ in L > 3 and at least 2 volumes 

■ m™ m > 400MeV ■ otherwise, or if (L min < 2fm) 

Moreover, the FLAG rules stipulate that results which are classified with at least one red tag 
and/or without a journal reference be excluded from global estimates. As an example we show the 
compilation of results for the strange quark mass from the FLAG report in Fig. 1. One observes 
that lattice estimates for m s obtained with Nf = 2 + 1 flavours of dynamical quarks appear to be 
somewhat smaller compared to the two-flavour theory, although this may be attributed to the fact 
that the results for Nf = 2 are typically older and may be more strongly affected by other systematic 
effects. A striking feature of the plot is that lattice estimates are broadly consistent, despite the fact 
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Table 1: Results for the strange quark mass, m^ s (2GeV), and the ratio m s /m uc j from simulations near the 
physical pion mass. Results marked by a dagger were obtained using mass reweighting. Numbers in red are 
based on non-perturbatively determined renormalisation factors. 

that they have been obtained for several different discretisations of the quark action. Moreover, 
the quoted uncertainties are much smaller than those attributed to sum rule results and the PDG 
average PDG2010. The FLAG report quotes the following global estimates, based on the results of 
refs. [9,10] and [11]: 

m^(2GeV) = 3.43(11) MeV, mf S (2GeV) = 94(3)MeV, m s /m ud = 27 .4 ±0.4. (2.1) 

This high level of accuracy raises the question whether systematic effects, in particularly those 
associated with the chiral extrapolation, are indeed controlled. 

We begin by discussing two recent calculations which do not rely on chiral extrapolations. The 
PACS-CS Collaboration [8, 12, 13] has used non-perturbatively 0(a) improved Wilson fermions 
and the Iwasaki gauge action at a fixed value of the lattice spacing to determine the light quark 
masses at the physical pion mass via mass reweighting. The BMW Collaboration [7, 15] has 
performed simulations with smeared tree-level improved Wilson quarks at pion masses as low as 
120 MeV. Quark masses were obtained by an interpolation to the physical pion mass. Results and 
some simulation details are shown in Table 1 . 

At m™ m = 156 MeV the PACS-CS Collaboration is almost at the physical point. Via a short 
chiral extrapolation, PACS-CS have obtained the results shown in the first row of Table 1 [12]. In 
a subsequent work they have proceeded to simulate with hopping parameters Kr*) ex t which, 
according to the chiral extrapolation of the results in [12], correspond to the physical pion mass. 
In order to compensate the slight observed mismatch between the targeted and actually measured 
pion mass, they have reweighted their ensembles using the single-histogram method [14]. The sec- 
ond line in Table 1 indicates that the resulting estimate for the strange quark mass is quite different 
from the value obtained from an extrapolation. However, this increase can be largely attributed 
to the use of non-perturbative renormalisation factors in [8] which were found to be 30% larger 
than their perturbative counterparts in ref. [12]. This is consistent with the observation that the 
difference between extrapolation and reweighting is much less pronounced for the ratio m s /m uc j 
in which the renormalisation factors cancel. One concludes that reweighting allows one to avoid 
chiral extrapolations at the expense of incurring a larger statistical error. Despite the large over- 
all uncertainty, the reweighted results by PACS-CS do not agree well with the global estimates 
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Figure 2: The effect of cuts applied to the lower limit of the pion mass interval on the light quark mass 
(left) and the ratio m s /m ul [ (right) [16]. Red error bars denote the systematic uncertainty. The vertical bands 
represent the global results of eq. (2.1). 



ineq. (2.1). This may be explained by the presence of other systematic errors, most notably lattice 
artefacts and finite-volume effects, which are not yet sufficiently controlled. 

For Wilson-type fermions, the BMW Collaboration has realised the lowest pion masses so 
far, at least at the three coarsest lattice spacings. However, their results for the light quark masses 
were not included in the global estimates according to the FLAG rules, since [7, 15] had not been 
published by the time when ref. [1] was completed. BMW have performed a comprehensive in- 
vestigation of systematic effects, including studies of the stability of their results under variations 
in the ansatz for the fit function. To this end they compared chiral fits based on either SU(2) ChPT 
at NLO or on a Taylor expansion, viz. 

m " d = 2B{ ~ 2 (47if]2 ln A|j ( 1+c - sA ) versus m ud = Ci+c 2 ml + C3ml + c 4 A, (2.2) 

where A parameterises the deviation from the physical strange quark mass. BMW have also applied 
cuts on the pion mass, by limiting its maximum value to either 340 or 380 MeV. 

Given the availability of ensembles with pion masses as low as 120 MeV, it is interesting 
to study the question of the impact of interpolations in the pion mass, compared to relying on 
extrapolations. In other words, how do the results for quark masses and their uncertainties vary 
when the lower limit of the pion mass interval to which the chiral fit is applied, increases gradually? 
Figure 2 shows the results for the light quark mass m uc / and the ratio m s /m uc i obtained after applying 
mass cuts at m % = 120, 200 and 240 MeV, respectively [16]. One clearly sees that the results 
are quite stable and consistent within errors, indicating that chiral extrapolations are under good 
control. Unsurprisingly, the error increases for longer extrapolations. A look at Table 2 shows that 
the error budget is, in fact, increasingly dominated by the fit ansatz when the minimum pion mass 
is shifted to larger values. It would be helpful to discount the possibility that these findings are 
obscured by lattice artefacts, since below-physical pion masses have only been simulated for the 
three coarsest lattice spacings. As a suggestion for an improved future analysis of the BMW data, 
the effects of imposing lower mass cuts should be investigated after the data have been extrapolated 
to the continuum limit for fixed values of m % . The above discussion shows that simulations at or 
below the physical pion mass allow for a systematic investigation into the quality of chiral fits. In 
order to rightfully claim an overall accuracy of a few percent in lattice estimates of the light quark 
masses, minimum pion masses of 250 MeV appear to be sufficient. 
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Table 2: Error budget for the light quark mass rri^f (2 Ge V) after applying cuts to the minimum pion mass of 
120, 200 and 240 MeV. The first row corresponds to the original results in [7]. The last six colums represent 
the relative contributions of individual systematic effects to the overall systematic error <7 syst . 



3. Systematics of chiral fits: meson decay constants 

Masses and decays constants of pseudoscalar mesons belong to the set of quantities whose 
dependence on the quark mass has been studied most extensively. Chiral fits using lattice data and 
the expressions of ChPT give access to the effective coupling constants (low-energy constants - 
LECs) of ChPT. These include the pion decay constant in the chiral limit, /, the quark condensate 
E, and also some of the LECs which enter at NLO in the chiral expansion (e.g. £3 and £4). Fur- 
thermore, the decay constants of the physical pion and kaon, f K , fx, as well as the ratio fx/ fx w& 
determined in this way. While the individual decay constants are often used to set the lattice scale, 
the ratio fa/fx is important for constraining the ratio | V t( .v|/|Vud| of CKM matrix elements. Both 
aspects will be covered in this section. The discussion here is restricted to lattice calculations in 
the /^-regime. 

Lattice results for the ratio fx./ fx are i n general quite stable and consistent among different 
collaborations. Examples of chiral extrapolations are shown in Fig. 3. The FLAG report provides 
separate global estimates for QCD with Nf = 2 and 2 + 1 flavours, i.e. 

/k//» = 1.193 ±0.005 (N f = 2+l), fx/ fit = 1.210 ± 0.006 ±0.017 (N f = 2). (3.1) 

The result for Nf = 2 + 1 is based on refs. [18, 20, 21], while the two-flavour result is identical to 
the value quoted in [22]. 

As was pointed out by Marciano [23], a precise lattice estimate of fx/ fx in conjunction with 
accurate experimental measurements of the leptonic decay rates of K^ 2 decays provides a stringent 
constraint on the ratio |V„i|/|Kt/| of CKM matrix elements, since 

TiK^nv^y)) \Vus\ 2 fjmx 
T{%^eUr)) ~ \Vud\ 2 ftm n K ' ' 

An additional constraint on \V m \ is provided by the form factor f+(q 2 ) which appears in the expres- 
sion for the rate of the decay K — > nlv . Lattice calculations for f+(q 2 ) are, consistent and equally 
precise compared with the effective field theory result of ref . [24] in which / + (0) was determined 
by invoking the Ademollo-Gatto theorem [25], which states that the corrections due to isospin and 
SU(3) flavour breaking are second order. Global lattice estimates for /+(0) are quoted in [1] as 

/ + (0) = 0.9597 ±0.0038 (N f = 2+l), / + (0) = 0.9604 ±0.0075 (Nf = 2) . (3.3) 

The current accuracy of lattice results for fx/ fx and /+(0) allows for a precision test of first-row 
unitarity of the CKM matrix, i.e. 

\V ud \ 2 + \V us \ 2 + \V ub \ 2 = l. (3.4) 
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Figure 3: The ratio /k/A plotted versus the squared pion mass from refs. [10, 12, 18, 19]. Open circles 
denote the extrapolated values at the physical pion mass which are shifted for clarity. The horizontal lines 
represent the global estimate from the FLAG report [1] for Nf — 2 + 1. 

The contribution from the ^-quark can be dropped, since \V uh \ 2 = 0(l(r 5 ), which is below the 
relevant level of accuracy in the following discussion. The experimental branching fractions are 

V us f K 



\V us \f+(0)= 0.2163(5) 



Vudfn 



0.2758(5). 



(3.5) 



Combining these values with the global lattice estimates for /Vf = 2 + 1 in eqs. (3.1) and (3.3) yields 



|^| 2 + |K.v| 2 = 1-002 ±0.015. 



(3.6) 



The precision of this test can be considerably enhanced by including another experimental con- 
straint, namely the determination of V ut \ from super-allowed nuclear j3-decays. Using the lattice 
result for the ratio /k/ fn which fixes | V„. s |/|V t( d| gives [1] 



\Vud\ 



2 + \V I 2 

1 I * us I 



0.9999 ±0.0006. 



(3.7) 



In this way, first-row unitarity is confirmed with per-mille accuracy, using experimental information 
and lattice estimates alone. The unitarity test is equally precise if the lattice estimate is provided 
by / + (0) instead of fa/fit- 

The above discussion suggests that lattice calculations of pseudoscalar meson decay constants 
are under very good control. Table 3 contains a compilation of recent estimates for decay con- 
stants, the hadronic radii ro and r\ and certain combinations thereof. Although the ratio /k//w is 
consistent among different calculations within the quoted errors, this is not necessarily true for the 
absolute values of decay constants and hadronic radii. On the assumption that there are no sig- 
nificant differences between two- and three-flavour QCD within the presently quoted errors, one 
finds that the ro determination from ETMC [22] (which uses the physical value of f n to set the 
scale) contradicts the estimate quoted by RBC/UKQCD [10]. On the other hand, RBC/UKQCD, 
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Table 3: Results for decay constants and the hadronic radii ro and r\ from simulations with Nf = 2 + \ 
dynamical flavours (RBC/UKQCD, PACS-CS, MILC) and N f = 2 (ETMC). The renormalisation of the axial 
current in [12,26] was done perturbatively. 
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Figure 4: Chiral extrapolations of f K from RBC/UKQCD. The left panel shows the comparison of SU(2) 
ChPT with an analytic ansatz, using a minimum pion mass of 290MeVat a w O.llfm [10]. In the right 
panel two additional points with m K = 170 and 250 MeV at a w 0.014fm have been included [27]. 

who determine the lattice scale from the mass of the Q. baryon, find a value for f % which is smaller 
than the experimental value. Similar observations apply to the results for /k, n and /k7i quoted 
by RBC/UKQCD and MILC [26]. 

One may suspect that the observed differences in the absolute values of /k and f n are linked to 
the normalisation of the axial current. In fact, only the correctly normalised matrix elements can be 
expected to approach the continuum limit with a rate proportional to the leading lattice artefacts. 
For quantities such as ro and r\ the situation is not much better, because little is known about 
the chiral behaviour one is to expect. In view of the importance of decay constants and hadronic 
radii for the overall scale setting, one must make an effort to understand the observed differences. 
Several collaborations have reported new results for these quantities at this conference [27-32]. 

RBC/UKQCD have investigated why the individual decay constants f K and /k are lower than 
experiment, while their ratio agrees with other simulations. To this end they have supplemented 
their existing data sets by two more ensembles with pion masses of 250 and 170 MeV [27]. In or- 
der to keep m™ m L > 4 for L/a = 32, lower pion masses could be simulated at the expense of having 
to use coarser lattice spacings (a « 0.14fm) compared to [10]. Figure 4 shows the impact of the 
additional data points on the chiral extrapolation of f n . Clearly, the ambiguity associated with the 
ansatz for the chiral behaviour is reduced: Extrapolations based either on SU(2) ChPT or on a Tay- 
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Figure 5: Chiral extrapolation of /k by ALPHA [31]. Left: schematic view of the paths of two different 
extrapolations in the (m„d,m2)-plane. Right: chiral extrapolations of o/k/Za in the variable y\ = m^/8n 2 f^ 
for three different bare couplings. 

lor expansion produce results for f n which agree very well within errors. However, while the ratio 
fit/ fx is consistent with the earlier result, the new and preliminary value of f % = 125(2)(3)MeV 
still appears to be smaller than the experimental value. It should be kept in mind that the entire 
range of pion masses which RBC/UKQCD have access to, involves different lattice spacings. The 
systematics of chiral fits can only be investigated reliably provided that the dependence on the 
lattice spacing is well understood. 

The consistency of different chiral fit ansdtze has also been studied by the ALPHA Collabora- 
tion for O(a) improved Wilson quarks [31]. In order to check the robustness of the extrapolation of 
fx to the physical pion mass, they have compared two different fit strategies. Denoting the hopping 
parameters by K\ and Kz, where JCi = K" sea , the first strategy amounts to adjusting K% until nv\/ f^ 
is equal to the experimental value. Repeating this for each sea quark mass defines, at leading or- 
der, a sequence of data points with m s + m u d = const., which can be extrapolated to the physical 
pion mass using the expressions from partially quenched ChPT In the second strategy the strange 
quark mass is held fixed: at each value of K" sea , the hopping parameter K2 is tuned such that the 
PCAC mass reproduces the fixed value of /I = m s . The resulting values of /k can then be extrapo- 
lated in the pion mass using SU(2) ChPT at NLO. A schematic view how the physical point in the 
(m uc i,m s ) -plane is approached with the two strategies is shown in Fig. 5 (left). The right panel of 
the figure compares the chiral extrapolation of the (bare) /k for three different values of j3. Even 
though they are based on quite different variants of ChPT, both strategies converge to the same 
values at the physical pion mass, which enhances the credibility of the chiral fits. After taking the 
renormalisation of the axial current into account the results can be used to set the lattice scale using 
/k. In this way ALPHA find that the three j8 -values correspond to a ^ 0.075,0.066 and 0.049 fm, 
respectively. These preliminary results agree well within errors with the scale determination via 
the mass of the Q. baryon [35] on the same ensembles. 

Decay constants are often used to calibrate the hadronic radii ro and r\, by determining com- 
binations such as f K ro at the physical pion mass in the continuum limit. Chiral extrapolations of ro 
and r\ are usually based on the assumption that they can be described by a polynomial in ra|. A 
detailed investigation of the chiral behaviour of ro was presented at this conference [32] (see also 
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Figure 6: Left: pion mass dependence of ro/ ro| ref computed on the CLS ensembles [32] (red symbols) 
compared to the results from ETMC [34] (blue symbols). The quantity ro | re f denotes the value at a reference 
pion mass. Right: continuum extrapolation of the slope parameter s in eq. (3.8). 

ref. [33]). In Fig. 6 the pion mass dependence of ro computed on the CLS ensembles with 0(a) 
improved Wilson quarks is compared to the results by ETMC [34]. As can be seen from the left 
panel of the figure, ro indeed shows a linear dependence on the squared pion mass at a given value 
of the bare coupling. By performing a linear fit according to 

ro/ r olref = A + sx («Vo) 2 , (3.8) 

one can determine the slope parameter s, which is found to depend quite strongly on the lattice 
spacing in the case of ETMC. The right panel of Fig 6 shows the extrapolation of s to a = 0. The rate 
with which this quantity approaches the continuum limit is indicative of the size of lattice artefacts 
of 0(m q a ). While the slope s from the two collaborations agrees very well in the continuum limit, 
there are sizeable corrections proportional to m q a 2 for twisted-mass QCD. This discussion serves 
as a reminder that the problem of mass-dependent lattice artefacts, which are formally of 0(a 2 ) 
must be addressed, since they will affect the determination of f n r$ and hence the calibration of ro. 
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4. Nucleon axial charge 

While hadronic uncertainties in the meson sector could be brought well under control, the 
situation for baryonic quantities is much less satisfactory. Despite many years of dedicated ef- 
fort, lattice results for nucleon form factors or moments of structure functions fail to reproduce the 
experimental values within the quoted uncertainties [36,37]. A prominent example is the axial 
charge, gA, of the nucleon. Lattice simulations using pion masses m K > 250 MeV typically under- 
estimate this quantity by 10 — 15 %. Even more worrisome is the observation that the gap is stable, 
i.e. the data show little if no tendency to approach the physical value as the pion mass is decreased. 
There is a broad concensus that uncontrolled systematic effects must be held responsible. 

The axial charge is an ideal observable to study lattice systematics for baryonic quantities: 
It is defined in terms of a transition matrix element at zero momentum transfer and hence the 
underlying kinematics is very simple. Second, it can be determined without the evaluation of 
quark-disconnected diagrams. Among the common sources of systematic error are lattice artefacts, 
the related issue of the correct normalisation of the axial current, and the influence of finite- volume 
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effects, which are known to be larger for baryonic systems. An obvious question is whether the 
chiral behaviour is sufficiently controlled in the calculations performed so far, or whether much 
smaller pion masses are required in order to make contact with the experimental value. Another 
issue which has received quite some attention recently, is the possible contamination of baryonic 
three -point correlation functions by contributions from higher excited states. This seems plausible, 
since the noise-to-signal ratio in baryonic correlation functions is much worse than for mesons. 
Thus, one cannot firmly rule out the possibility that excited state contributions are still present 
within the relatively short Euclidean time interval before the signal is lost. 

The theoretical foundations of baryonic ChPT, which is used to constrain the chiral behaviour 
of gA, are unfortunately on a weaker footing compared to the mesonic sector. Since the mass 
gap between the nucleon and the nearest resonance, i.e. the A, is much smaller than the mass 
scale defined by the nucleon itself, it is difficult to define a consistent chiral counting scheme. 
The established formalisms include Heavy Baryon ChPT [38], the infrared regularisation of loop 
integrals [39] and the related extended on- mass shell regularisation [40], some of which have been 
carried to high orders in the expansion. Another approach is the so-called small-scale expansion 
(SSE) [41], in which the nucleon-A splitting is treated as a small parameter and included in the 
chiral power counting in the framework of Heavy Baryon ChPT. One severe drawback for the 
interpretation of lattice data is the large number of coupling terms, each of which carries a low- 
energy constant. Some of these LECs can be constrained from phenomenology, but unless one has 
access to extremely detailed information from lattice simulations deeply in the chiral regime, it 
seems impossible to determine the full set. 

The bare value of gA can be extracted from a suitable ratio of two- and three-point functions. 
In the simplest case, i.e. when the same operators are used to create and annihilate the nucleon, the 
expression reads 

R A (t,t s ) = f ' ( '^ >>0 g b r + 0(e-^) + 0(e- A ^). (4.1) 

Here t s denotes the Euclidean time separation between the initial and final nucleons, while the 
axial current is inserted at time t with < t < t s . Due to the rapidly increasing statistical noise, 
typical values of t s are of order 1 fm, and thus the correlation functions must reach their asymp- 
totic behaviour for separations t, (? s . — t) < 0.5 fm. Since the gap An between the nucleon and its 
first excitation is expected to scale like An ~ 1m % in the chiral regime, it is clear fromeq. (4.1) 
that corrections from excited states become increasingly important as the physical pion mass is 
approached. 

Several collaborations have investigated the issue of excited state contaminations recently. 
Using the CLS configurations generated with Nf = 2 flavours of 0(a) -improved Wilson fermions, 
the Mainz group has calculated baryonic three-point functions for several different source-sink 
separations t s [42]. After computing the so-called "summed insertions" [43] according to 

S A (t s ) = Y,R A (t,t s ) ^ const. +gTt s + 0(t s e-*» t '), (4.2) 

they determine the axial charge from the linear slope of Sa in the source-sink separation t s . Since 
t s > t, (? s — t) by construction, it is clear that the corrections due to excited state contamination in 
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Figure 7: Left: preliminary results for gA computed on the CLS configurations using the technique of 
summed insertions [42]. Right: compilation of recent results for gA from refs. [42,46,47]. The experimental 
result for gA is represented by the pink asterisk. 



eq. (4.2) are parametrically more strongly suppressed than for the simple ratio R&(t,t s ). Other ways 
to address excited state contamination include the use of multi-exponential fits [44] and systematic 
studies of the dependence of gA on the source-sink separation for t s as large as 1.9fm [45]. All 
these efforts do not allow for a firm conclusion at this stage. The preliminary results by the Mainz 
group (left panel of Fig. 7) suggest that summed insertions lead to a better agreement with the 
experimental result for g\. ETMC [45] report the absence of a bias in gA at m K = 380 MeV but 
see some evidence for a distortion in the case of (x) u -d- For the latter quantity, LHPC [44] can 
confirm that multi-exponential fits lead to a better agreement with experiment as the pion mass is 
lowered, albeit with a larger statistical error. 

New results on gA by the RBC/UKQCD Collaborations were presented at this conference 
[46], computed on the set of gauge configurations which included the recently added lighter pion 
masses discussed in section 3. They report a stable gap between their preliminary results and 
the experimental value of gA across the entire mass range. The favoured explanation offered by 
RBC/UKQCD is that the discrepancy is a result of finite-volume effects rather than excited state 
contamination. A compilation of recent results for gA [42, 46, 47] is shown in the right panel 
of Fig. 7, where the chiral behaviour is compared among different groups after applying the cut 
m K L > 4. Despite the very different systematics concerning the discretisation of the quark action, 
the values of the lattice spacing and the numerical procedures to extract gA from the measured cor- 
relation functions, the results are broadly consistent with each other. However, it would appear that 
those calculations which address the issue of excited state contamination compare more favourably 
with the experimental value. 



5. Summary and conclusions 



Chiral extrapolations have been a persistent source of systematic errors in lattice calculations. 
In this review I have tried to assess the reliability of chiral extrapolations in order to investigate 
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whether the claimed accuracy of lattice results for several phenomenologically interesting quan- 
tities is justified. The emergence of simulation data around the physical pion mass was crucial, 
since it allowed for a systematic study into the effects of replacing the chiral extrapolation by an 
interpolation. 

Pion masses in the range of 250 — 400 Me V appear to be sufficient to guarantee that lattice 
estimates for the light quark masses can be obtained with overall uncertainties at the level of a few 
percent. Similarly, the chiral behaviour of the ratio /k/A i s under good control. The latter allows 
for a precise determination of the ratio \V US /V u d\ and for a test of first-row unitarity with permille 
accuracy, based on lattice results and experiment alone. For individual decay constants, however, 
small inconsistencies among different calculations remain and must be resolved. The separation of 
lattice artefacts from systematic effects associated with the description of the chiral behaviour must 
be improved not only for /k and f n but also for quantities such as r§. 

In spite of these successes, one finds that lattice calculations for the axial charge are still in an 
unsatisfactory state, since the chiral behaviour of gA is clearly obscured by systematic effects. With 
the presently available data it is difficult to decide whether one single cause is chiefly responsible or 
whether it is a convolution of finite-volume effects, excited state contamination and lattice artefacts. 
Due to the unfavourable signal-to-noise ratio of baryonic correlation functions it is likely that this 
can only be resolved via an enormous increase in statistics. 
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